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Abstract: Using a renormalized linked-cluster-expansion method, we have extended to order 
f} 23 the high-temperature series for the susceptibility \ an d the second-moment correlation length 
£ of the spin-1/2 Ising models on the sc and the bcc lattices. A study of these expansions yields 
updated direct estimates of universal parameters, such as exponents and amplitude ratios, which 
characterize the critical behavior of \ an d £• Our best estimates for the inverse critical temperatures 
are (3 s a c = 0.221654(1) and /3 b c cc = 0.1573725(6). For the susceptibility exponent we get 7 = 1.2375(6) 
and for the correlation length exponent v — 0.6302(4). The ratio of the critical amplitudes of \ 
above and below the critical temperature is estimated to be C+/C- = 4.762(8). The analogous 
^ . ratio for £ is estimated to be /+//- = 1.963(8). For the correction-to-scaling amplitude ratio we 

O ' obtain af /at = 0.87(6). 

PACS numbers: 05.50+q, 11.15.Ha, 64.60. Cn, 75.10.Hk 
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I. INTRODUCTION 
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As a part, of an ongoing long-term program of computer-based calculations and analyses of HT series for two- 
dimensionalu and for three-dimensionarm lattice spin models, we have extended by two terms the high-temperature 
(HT) series for the spin-1/2 Ising model on the simple-cubic (sc) and the body-centered-cubic (bcc) lattices. In the 
first analysis presented here, we shall restrict to consider the HT expansions through /3 23 for the susceptibility x an d 
the second-moment correlation length £, mainly in order to update the direct estimates of the corresponding critical 
0$ ■ indices 7 and v. 

For the sc lattice, the longest^xpansions of these quantities already in the literature reach order /3 21 _.They were 



Qh obtained and analyzed in Refsrm only a few years ago. In the case of the bcc lattice, the published serie sMfor 

£, also extending through /3 21 , were calculated by B.G. Nickel two decades ago. The progress in such computations has 
been slow due to the exponential growth of their complexity with the order of the expansion_so that even adding only 
a few terms to the present results is a laborious task. Within the renormalized linked-clusteiEl expansion method, used 
in our work, one must overcome many problems of combinatorial nature concerning graph generation, classification 
and partial resummation, and a special effort must be devoted to keep under strict control the numerous possible 
sources of error. In our case, a final severe test is provided by having the program to reproduce, in three dimensions- 
established data like the series for the nearest-neighbor spin correlation on the sc lattice, which is already|-t,abulatedLl 
through /3 27 , and_in two dimensions, the series for \ and [ii on the. simple square lattice, which are knowntffl through 
/3 35 and beyondfltl After the completion of this work, a preprintlil] has been issued which also reports independently 
extended expansions for \ an d on both-the the sc and the bcc lattice through orders /3 23 and /3 25 respectively. Our 
series coefficients agree with those of Ref.EHl as far as the expansions overlap. This adds further confidence about the 
correctness of the results since our implementation of the linked cluster expansion procedure is rather different from 
that described in Ref.El 

Any enrichment of the exact informations on the 3d Ising model is still of general interest. Here we have used these 
novel data to improve the knowledge of the nonnegligible singular correctionslliJ to the leading critical singularities of 
X and ^^nd, as a consequence, the accuracy of the direct HT series estimates of all critical parameters. As stressed 
in Refrrtil, the corrections to scaling first showed up unambiguously when the bcc series were extended to order /3 21 , 
the last three coefficients being crucialcl. It is therefore helpful to produce more coefficients, in order to stabilize and 
possibly refine the quality of the information extracted from the series. 

The plan of this note is as follows: after setting our notational conventions in Sec. II, we tabulate the series 
coefficients for \ and ^ 2 through order 23, with respect to the usual HT expansion variable v — th(fl). In Sec. 
Ill we report the results of our extrapolations for the critical temperatures, for the critical exponents 7 and v, for 
the universal ratio C+/C_ of the critical amplitudes of the susceptibility above and below the critical point, for 
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the analogous ratio /+ / /_ of the correlation- length amplitudes and for the ratio /a+ of the correction-to-scaling 
amplitudes^. Our estimates are compared with the latest numerical calculations by series, by stochastic methods 
and by perturbative teppcmalization group (RG) techniques, in the fixed-dimension (F-Di) approacht3~0 and in the 
e-expansion approacHi!rE3. Less recent studies have been already reviewed in our Refsrm. 

II. DEFINITIONS AND NOTATIONS 

In order to introduce our notation, we shall specify by the Hamiltonian 

{x,x') 

the nearest-neighbor three-dimensional spin-1/2 Ising model in zero magnetic field. 

Here s(x) = ±1 is the spin variable at the lattice site x, and the sum extends over all nearest neighbor pairs of 
sites. We shall consider expansions in the usual HT variable j3 = J /ksT called "inverse temperature" for brevity. 
However, for convenience, we shall tabulate the series coefficients with respect to the expansion variable v = th(/3). 

The susceptibility is expressed in terms of the connected two-spin correlation function (s(x)s(y)) c by 

oo 

#) = ^( S (0) S (f)) c = l + ^a#; (2) 

x r=l 

and the second moment of the correlation function is defined as 

oo 

M =^f 2 ( S (0) S (f)) c = ^6 r /3'-. (3) 

x r=l 

In terms of \ an d 1^2 the second-moment correlation length £ is defined by 

em = l§ y w 

For easy reference we report here the complete expansions of x an d (J-2, rather than only the lastly computed two 
coefficients. For the susceptibility on the sc lattice we have 

X sc (v) = 1 + 6v + 30v 2 + 150v 3 + 726v 4 + 3510u 5 + 16710i> 6 + 79494i; 7 + 375174u 8 + 1769686u 9 + 8306862u 10 
+38975286D 11 + 182265822t> 12 + 852063558ti 13 + 3973784886u 14 + 18527532310t> 15 + 86228667894i/ 6 



4-864404776466406-u + 4006394107568934tj + . . . 

for the second moment on the sc lattice: 

fi^iv) =611 + 72ti 2 + 582-u 3 + 4032-u 4 + 25542t> 5 + 153000« 6 + 880422t> 7 + 4920576-u 8 + 26879670u 9 + 144230088« 10 
+762587910U 11 + 3983525952t> 12 + 20595680694i> 13 + 105558845736« 14 + 536926539990t> 15 + 2713148048256« 16 



for the susceptibility on the bec lattice: 

X bcc (v) = 1 + 8v + 56v 2 + 392v 3 + 2648u 4 + 17864i> 5 + 118760i> 6 + 789032u 7 + 5201048t> 8 + 34268104u 9 + 224679864U 1 
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+ 1472595144D 11 + 9619740648u 12 + 62823141192w 13 + 409297617672u 14 + 2665987056200u 15 + 17333875251192-u 16 
+ 112680746646856« 17 + 731466943653464u 18 + 4747546469665832t> 19 + 30779106675700312u 20 + 199518218638233896t> 21 
+1292141318087690824« 22 + 8367300424426139624i> 23 + . . . 

for the second moment on the bcc lattice: 

/4° c («) = 8v + 128v 2 + 1416d 3 + 13568t> 4 + 119240t> 5 + 992768t> 6 + 7948840u 7 + 61865216« 8 + 470875848i> 9 + 35219548161) 11 
+25965652936D 11 + 189180221184ti 12 + 1364489291848u 13 + 9757802417152u 14 + 69262083278152-u 15 + 488463065172736u 16 
+3425131086090312-U 17 + 23896020585393152u 18 + 165958239005454632-u 19 + 1147904794262960384« 20 
+ 7910579661767454248u 21 + 54332551216709931904« 22 + 372033905161237212392?; 23 + . . . 



III. ANALYSIS OF THE SERIES 



In terms of thejeduced inverse temperature r# = 1 — (3/ p§ , the asymptotic critical behavior of the susceptibility 
is expected to bellil 

X *{(3) c C*(T*)-l(l + a+#(r#) e + ... + e+*r# + . . . ) (5) 

as the critical point (3f is approached from below. (Here and in what follows, the superscript # stands for either 
sc or bcc, as appropriate, and will be dropped whenever unnecessary. The index +(-) denotes, as usual, quantities 
associated with the high (low) temperature side of the critical point.) Similarly, for the correlation length £, we expect 

£#(/?) ~ f# {t*)-» (l + 4*(r*f + ... + e + *r* + . . . ) (6) 

as r — > + . 

The exponents 7, v and 9 are universal quantities, whereas the critical amplitudes C*, the amplitudes 

of the leading nonanalytic correction-to-scaling terms and the amplitudes e+*, ef of the leading analytic 
corrections are nonuniversal, as suggested hv the superscript #. Experimentally accessible universal combinations 
can be formed out of the critical amplitudesEHl. Here we shall be concerned with series estimates of the universal ratios 
C + /C_, /+//- and a^/a+. Notice that for the critical amplitudes we have adopted the notation of Ref.ta and of 
other recent studies rather than that of Ref.El. 



A. Estimates of the critical points 



As a first step of the analysis, we shall examine the series for the susceptibility whose coefficients have the smoothest 
pattern of behavior, so that they are generally used to estimate the critical temperatures. These estimates will also be 
used to bias the determination of the critical exponents and of the universal amplitude ratios, therefore their accuracy 
is crucial. Let us .begin by considering the results obtained by a very efficient variant of the ratio method introduced 
by J. Zinn-JustinEij, (see alsda). 

We evaluate [3 C from the sequence 

ou. = (^^/w S : +Sn - 2 j = Pc + o(4^) (7) 

a n a n -i 2s n (s n -s n -2) n 1+a 



where 



(in(A^)-i +i n ( a "-3 yi\ /2 ( 8 ) 
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This is an unbiased method, in the sense that no additional accurate information must be used together with the 
series in order to get the estimates of the critical parameters, but we found useful to improve the procedure by biasing 
it with the value of 9 as follows. For sufficiently large n, the sequence of estimates {f3 c ) n shows very small regular 
oscillations due to the loose structure of the lattice. Moreover the odd and even subsequecences of (f3 c ) n have a 
residual decreasing trend which is very nearly linear on a l/n 1+0 plot, as suggested by eq. (Q). Therefore, simply 
taking the highest order term of the sequence ((3 c )n as the final estimate, would be an inadequate choice. We have 
preferred to extrapolate separately to n — > oo the successive odd and even pairs of estimates (/3 c )n, assuming that we 
know the value of 9 well enough. The two sequences of extrapolated values need further extrapolation which allows 
also for the small residual curvature of the plot and leads to the final estimates j3 s c c = 0.221654(1), in the case of the sc 
lattice, and = 0.1573725(6), in the case of the bcc lattice. The errors we have reported, account generously both 
for the present uncertainty in 9 (whose effects in this analysis are very small anyway) and for the uncertainty of the 
second extrapolation—For the correction-to-scaling exponent we have assumed the value 9 — 0.504(8), obtained by the 
FD perturbative RGO. Also in the rest of this note the central values of all ^-biased estimates will refer to this value. 
However, in the calculations of this and the next subsection, we have also considered a much larger uncertainty in 
order to make sure that oup-^esults are compatible with somewhat higher central values su&b,as 9 = 0.52(3), proposed 
in Ref.B (as well as in Ref.cil, with a smaller error), or with 9 = 0.53(1) suggested in Refri, An even larger central 
value 9 — 0.54(3) was indicated in Refs.E3'Ell, while an experimental measure reported in Ref.E2l yields 9 — 0.57(9). In 
the case of the bcc lattice, as an example of our extrapolation procedure, we have reported in Table 1 the last eight 
terms of the sequence (f3 c ) n & n d the results of the initial extrapolation of the last six successive alternate pairs of 
terms. Our final result for the critical inverse temperature of the Ising model on the sc lattice is completely compatible, 
although much less precise than the value (3 S C C — 0.22165459(10) obtained from an extensive MonteCarlo (MC) study 
by a dedicated Cluster ProcessorcJ and generally considered as the best available estimate. Our central value of /3c CC , 
obtained similarly-is, paly slightly smaller, but more precise than the value-j3* cc = 0.157373(2) suggested in Nickel 
and Rehr analysisaSEil We should finally mention that Prof. D. StauffercJ kindly informed us that he still tends 
to favor the somewhat larger central estimate /3* c = 0.221659, basing on the HT .analysis in Ref.Eil, as well as on 
his own recent simulatior£3 of the critical dynamics and on analogous work in RefH3. We also recall that a similar 
value = 0.2216595(26) was indicated a decade ago in the MonteCarlo simulation of Ref.E3. In the context of our 
analysis, these values lie approximately halfway between the highest order approximant (/3;? c )23 ~ 0.221667 and our 
final estimate obtained from extrapolation. To close this section, three remarks are in order. First: the reliability 
of our analysis procedure has been corroborated, by repeating it with the recently computed (9(/3 26 ) series for the 
self-avoiding-walk (saw) model on the sc latticed. This is a relevant test because the structure of the corrections to 
scaling (namely the sign and size of the correction amplitude and the value of the confluent exponents) is expected 
to be quite similar to the Ising sc case. For the saw model we have observed that the central value for f3 c indicated 
by our procedure is essentially stabilized after reaching the order /3 23 and agrees closely with that indicated in Ref.O, 
while the error decreases as higher order coefficients are included in the analysis. Our procedure has also been tested 
and confirmed by other arguments in Ref.E-tl. Second: due to the higher coordination number of the bcc lattice, 
the corresponding series have a greater "effective length" than the sc series, and therefore all estimates obtained for 
the bcc lattice will be systematically more accurate. Third: as expected, the inclusion in our analysis of the two 
additional coefficients for the expansion of \ 011 the bcc lattice, computed in Ref.E2l, does not essentially modify our 
central estimate of (3^ cc , but only reduces its uncertainty to the value reported here. 



B. Estimates of the critical exponents 

By using a related variantEltd of the ratio method and by analogous arguments, fairly good estimates can be obtained 
also for the exponents 7 and v. We construct the approximation sequence 

2(s n + s n - 2 ) . 1 , , 

7 " = 1+ Gs„- s „_ 2 ) 2 =7 + 0( ^ ) (9) 

with the same definition as above for s n . Also in this case, for sufficiently large n, the successive estimates 7„, 
(as well as the analogous ones v n obtained from the series coefficients of £ 2 ), appear to be nearly linear on a l/n e 
plot, and therefore we can follow an extrapolation procedure completely analogous to the one previously described. 
However, in the exponent calculation, the corrections are a priori larger and therefore the procedure involves relative 
errors larger than in the case of (3 C . In order to illustrate this numerical procedure in the case of the bcc lattice, 
we have reported in Table 1 the last eight terms of the sequence j n and the results of the extrapolation of the last 
six successive alternate pairs of terms. The estimates inferred from the analysis of these data are: 7 = 1.2378(10) 
and v — 0.629(2) in the case of the sc series and 7 = 1.2373(6) v = 0.629(1) from the bcc series. As expected, the 
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relative uncertainties for the exponent v are larger because of the slower approach of the second moment series to its 
asymptotic behavior. The dependence of these estimates on the value of used in the extrapolation can be expressed 
as follows: 7 = 1.2378 + 0.016(0 - 0.504) ± 0.0010 and v = 0.629 + 0.02(6* - 0.504) ± 0.0020 in the case of the sc lattice; 
7 = 1.2373 + 0.012(61 - 0.504) ± 0.0006 and v = 0.629 + 0.016(6* - 0.504) ± 0.0010 in the case of the bcc lattice. 

In order to confirm these estimates for the exponents, we shall resort also to (unbiased and biased) analyses by 
inhomogeneous differential approximants (DA's)ElH. By unbiased DA's, we obtain somewhat larger estimates both 
for the critical inverse temperatures and for the exponents, which, however, show a clear decreasing trend. Therefore 
also these data should be further extrapolated, but, unfortunately, this is not as straightforward as in the case of the 
Zinn-Justin method. Thus we did not insist on this route and preferred to perform biased^eries analyses, either i) by 
the first-order simplified differential approximants (SDA) introduced and discussed in Ref.u, in which both (3 C and the 
correction-to-scaling exponent are fixed, or alternatively ii) by conventional second order inhomogeneous-DA's, in 
which 9 and (3 C are varied in a small neighbourhood of their expected values, following the method of Ref.tl Let us 
also add that in all cases in which we have relied on SDA's, we have,aiso repeated the same calculation, first subjecting 
the series to the biased variable change introduced by R. Roskiesca in order to regularize the leading correction to 
scaling and then computing simple Pade approximants. In this way we have always obtained completely consistent 
results, although they are sometimes affected by larger uncertainties. 

We have used the procedure i) to study the residue of the log-derivative of x ° r of £ 2 at the critical singularity. In 
the case of the sc lattice series, rather than our own estimate of (3 C , we have used the more accurate (but otherwise 
completely consistent) value (if = 0.22165459(10) of RefH Thus we estimate 7 = 1.2378(10) and v = 0.6306(8). 

In the analysis of the bcc lattice series, we have taken as a bias the value suggested by our extended ratio-method 
analysis (3 h c cc = 0.1573725(6). In this case we get the values 7 = 1.2375(6) and v = 0.6302(4). By using Fisher scaling 
lawEZl, we get rj = 0.037(3) from the sc series and rj = 0.036(2) from the bcc series. For both lattices we have used the 
same value (and uncertainty) of 9 as previously discussed and we have easily allowed for the residual decreasing trend 
of the exponent estimates, because SDA's values show a smaller spread than DA's. We can also mention that in the bcc 
lattice case, the linearized dependence of the exponent central estimates on the bias values of j3 c and 9 can be described 
as follows : 7 = 1.2375 + 0.01(0- 0.504) + 90.(/3 c - 0.1573725) and v = O.63O2 + O.O15(0-O.5O4)+4O.(/3 C -O.1573725). 

We shall take as our final estimates for the exponents those obtained by SDA's from the bcc lattice, which are best 
converged. r— ■ .— ■ 

The so-called M2 method of Ref.E3 is a very useful extension of the above mentioned Roskies' procedureE3. In the 
case of the bcc lattice it suggests j3^ cc = 0.1573720(4) with 7 = 1.2374(4) and 9 = 0.56(3), in good consistency with 
the other approaches. On the other hand, in the case of the sc lattice, the results of the M2 method at order /3 23 , 
namely-$;? c = 0.221659(2) 7 = 1.2395(5) with 9 = 0.50(2) are not essentially changed with respect to those obtained 
in Ref.Ei from the analysis of our previous 0{(3 21 ) series. 

In conclusion, provided that the sequences of estimates are carefully extrapolated using the independently computed 
value of 0, the determination of the exponents by the improved ratio method and by biased DA's or SDA's are 
completely consistent, though the latter method gives slightly more accurate results. At this order of expansion, 
asymptotic trends seem to be already stabilized and the uncertainties in the HT series estimates are significantly 
reduced. A sample of recent estimates of the critical exponents is reported in Table 2 and briefly commented in the 
rest of this subsection. The-agreement of our results with the values 7 = 1.2396(13) and v — 0.6304(13), indicated 
by the FD perturbative RGB, or the values 7 = 1.2380(50) and v = 0.6305(25), suggested by the e— expansiontH, is 
still good. However, we should observe that, in the years, as the length of the HT series has increased, the exponent 
estimates have been moving towards the slightly lower central values 7 « 1.237 and v w 0.630. indeed very similar- 
values had already been suggested some time ago by J.H. Chen, M. E. Fisher and B. G. NickelEj who studied soft 
spin models of the Ising universality class, chosen so to have negligible amplitudes for the leading corrections to 
scaling. Within this approach, the analysis^ of HT series through order /3 21 for the bcc lattice gave 7 = 1.237(2) .and 
v = 0.6300(15). A studyLJ of the HT series through 0(/3 20 ) for the sc lattice, along the same lines as in RefsB^, 
indicates 7 = 1.2371(4) and v = 0.63002(23). Recently, this methopLwas adapted also to MC simulations in Ref.E-3, 
which reports 7 = 1.2372(17) and v = 0.6303(6). Analogously in Ref.E3, the estimates v = 0.6296(7) and 77 = 0.0358(9) 
are obtained, implying 7 = 1.2367(20). Even lower central estimates of the exponents, namely 7 = 1.2353125) and 
v = 0.6294(10) have been obtained in a MC simulation of the Ising model by a finite-size scaling analysiso which 
allows for the corrections to scaling. 

C. Estimates of universal amplitude ratios 

By taking advantage also of the low-temperature expansion of x on the sc lattice, extendedjia RefEl to order u 26 , 
(here u = exp(— 4/3)), and of the older series for the bcc lattice computed to order it 23 in Refo, we can give a new 
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direct estimate of the universal ratio C+ / C_ . Using the low- temperature series for £ 2 , computed for the sc lattice in 
Ref.Cj through it 23 , we can also compute the ratio /+ C //- C - These quantities have been repeatedly evaluated in recent 
years by various techniques, with increasing accuracy. 

We have used first order SDA's to compute C± = lim r _o± l T l 7 X # - In the sc lattice case, by choosing f3^ c = 
0.22165459, 7 = 1.2375 and 9 = 0.5, we obtain Cp/CL = 4.762(8). Here the uncertainty refers to the sharp bias 
values above indicated. In order to compare this result with others obtained by slightly different assumptions, the 
dependence of our estimate on the bias values 76 and 9b for the critical and the correction exponents can be linearly 
approximated by C s ^/C s _ c = 4.762 + 6.06(1.2375 - 7 6 ) + 0.7(9 b - 0.5) ± 0.008. The ratio is insensitive to the choice of 
/3^ c within its quoted uncertainty. In the case of the bcc lattice, we have used the same bias values for 7 and 9 together 
with (3 bcc = 0.1573725, obtaining C+/C_ = 4.76(3). In this case the error (mainly coming from the uncertainty of 
C&_ cc ) j s l ar ger, but the result is completely consistent with the sc lattice estimate. 

The experimental measurements of this ratio range between 4.3 and 5.2E30E3 and are perfectly compatible with 
our estimates. Other recent numerical evaluations are summarized in Table 3. .Hjiwever some comments are helpful for 
understanding these results. The previous evaluations by M.E. Fisher and coILlZej, used shorter series and bias values 
7 = 1.2395 and (3 S C C = 0.221630 somewhat different from ours, thus yielding the slightly larger value 4.95(lJj 
ratio C+ / C_ can also be obtained feom approximate parametric representations of the scaling equation of state 
here we quote only the most repented such estimate: 4.77(2). 

The MC simulation of Ref.EJ gave the somewhat4arger value 5.18(35), while the more recent and higher precision 
study of RefJla yields 4.75(3) and the work of RefB reports 4.72(11). 

Within the &=expansion approach to the RG, the estimate 4.73(16) is obtained, while the FD expansion gives the 
result 4.79(10)0. The value 4.72(17) was obtained in RefM 

In a similar way, we have computed /+ c //- c = 1.963(8), assuming v = 0.6302. The quoted uncertainty allows also 
for the uncertainties in the estimates of v and 9. Other estimates appearing in [tip-recent literature are summarized 
in Table 4. Our result compares weLLwith the estimate 1.96(1) obtained in RefEO'Eil by shorter series as well as with 
the recent estimate 1.961(7) of RefM The MC estimate of Ref.N was 2.06(1), whereas in RefB the value 1.95/2) 
is reported. The latest e— expansion estimateta is 1.91 (with no indication of error bars) and the FD estimateca is 
2.013(28). The recent experimental estimates of this ratio range between 1.9(2) and 2.0(4). 

We believe that the close agreement between our series estimates and the latest determinations of these universal 
ratios is due to the careful allowance of the confluent corrections to scaling by SDA's. Indeed, even using the longer 
series presently available, simple Pade approximants, notoriously inadequate to describe the singular corrections to 
scaling, suggest estimates sizably larger, while the conventional DA's lead to a wider spread in the estimates. Further 
improvements of the direct series determination of these ratios should probably await for an extension of the low- 
temperature series. 

By using only the HT extend ed seri es presented here, we can also reevaluate the universal ratio a^~/a+. Let us 
recall that, as observed in RcfMclcaa and argued in earlier studied-^ for the spin-1/2 Ising model on the sc, the 
bcc and the fee lattices, the amplitudes of the leading correction-to-scaling terms have a negative sign, both for 
the susceptibility and the correlation length. The values of these amplitudes can be most simply determined, also 
in this case, by using the SDA's above mentioned. Our estimate for the universal ratio between these amplitudes: 
a^ sc /a+ sc — 0.95(15) from the sc lattice series, and a~^ bcc /a+ bcc = 0.87(6) from the bcc lattice series, improves the 
accuracy of our previous resultsi obtained from the analysis of shorter series. These results have to be.campared with 
the FD result a+/a+ = 0.65(5) obtained in Ref.E3 and with the HT result a+ bcc /a+ bcc w 0.85 of RefflB 

IV. CONCLUSIONS 

We have extended through order /3 23 the HT expansions of the susceptibility and of the second correlation moment 
for the spin-1/2 Ising model, on the sc and the bcc lattices. As a first application of our calculation, we have updated 
the direct HT estimates of universal critical parameters of the Ising model with some improvement over previous 
analyses in the accuracy and in the agreement with the latest calculations by approximate RG methods and by 
various numerical methods. 
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TABLE I. The sequences of approximants for j3 c and 7 defined by eq.(|7]) and eq.(|9|), respectively, and the sequences of the 
appropriate extrapolations using alternate pairs, as obtained from \ on the bcc lattice. For the extrapolations we have assumed 
that 6 = 0.504. 



_n (Pch from Eq.(pj ) Extrapol. of (/3 c ) n f 7l from Eq.(|g| ) Extrapol. of 7, 

18 0.1573815 1.244335 

19 0.1573806 1.244174 

20 0.1573807 0.1573761 1.244049 1.238519 

21 0.1573800 0.1573759 1.243889 1.238114 

22 0.1573799 0.1573743 1.243760 1.237595 

23 0.1573793 0.1573746 1.243620 1.237599 

24 0.1573791 0.1573739 1.243501 1.237475 

25 0.1573787 0.1573740 1.243374 1.237421 



TABLE II. A comparison among recent estimates of the critical exponents 7 and v. 



ci 
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Scries 
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M< 



M 



MC 



1.2375(6) 
0.6302(4) 



1.237(2) 
0.6300(15) 



1.2372(17) 
0.6303(6) 



1.2353(25) 
0.6294(10) 



FD-cxp 



1.2396(13) 
0.6304(13) 



1.2380(50) 
0.6305(25) 



1.2371(4) 
0.63002(23) 



1.2367(20) 
0.6296(7) 



TABLE III. A comparison among recent estimates of the susceptibility universal amplitude ratio C+/C- 



This work 
4.762(8) 



SerieOg 
4.95(15) 



Eq.Statotl' 

4.77(2) 



5.18(35) 
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M 



4.75(3) 
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